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The renormalization group approach in three dimensions is used to estimate the 
universal critical value of dimensionless sextic effective coupling constant for the 
Ising model. Four-loop expansion for g$ is calculated and resummed resulting in 
<7g = 1.596, while the most accurate estimate for g$ is argued to be equal to 1.61. 
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The critical thermodynamics of the three-dimensional Ising model is known to be de- 
scribed by Euclidean scalar field theory with the Hamiltonian 



H — d 3 x 



~ l -ml^ + + \^\ , (1) 

where a bare mass squared is proportional to T — Tjp\ being the phase transition 
temperature in the absence of the order parameter fluctuations. Taking fluctuations into 
account results in renormalizations of the mass mo — > m, the field cp = (p^^/Z, and the 
coupling constant A = mZ^Z^g^. Moreover, thermal fluctuations give rise to many-point 
correlations < ip(xi)tp(x2)...(p(x2k) > and, correspondingly, to higher-order terms in the 
expansion of the free energy ("effective action") in powers of magnetization M: 
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F(M) - F(0) = -m 2 Z- l M 2 + mZ^g^M 4 + £ m^ k Z~ k g 2k M 2k . (2) 

^ fc=3 

In the critical region, where fluctuations are so strong that they completely screen out the 
initial (bare) interaction, the behaviour of the system becomes universal and dimensionless 
effective couplings g 2 k approach their asymptotic limits g% k , i.e. they assume constant values 
which are also universal. 

Recently, several attempts have been made to find numerical values of higher-order 
coupling constants g$, g$, etc. at criticality; particular attention has been paid to the sextic 
effective interaction g 6 . C. Bagnuls and C. Bervillier solving exact renormalization group 
equations in the local potential approximation have obtained g$ = 2 AO M. M. Tsypin has 
performed extensive Monte Carlo simulations of the Ising model behaviour in an external 
magnetic field and found gl = 2.05 ± 0.15 ||. T. Reisz has deduced the estimate g$ = 
1.92 ± 0.24 analyzing linked cluster expansions for 0(n)-symmetric lattice model in the 
Ising case n — 1 0. 

These estimates obtained by means of essentially different methods are seen to be con- 
siderably scattered. Hence, the question about the numerical value of g§ remains, in fact, 



open. On the other hand, there is a theoretical tool which has proven to be very efficient 
when used to study the critical behaviour of various phase transition models. We mean 
the renormalization group (RG) approach in three dimensions. This method has enabled 
one to calculate critical exponents and fixed point locations with an accuracy of order of 1 
per cent for simple 0(n)-symmetric models [§]|)]|| as well as for more complicated systems 
possessing two and three fllOlJlll quartic coupling constants in their Landau- Wilson 

Hamiltonians. It is quite reasonable therefore to employ the field-theoretical RG approach 
in three dimensions for calculation of gl and other universal higher-order couplings. 

In this Letter, we estimate g^ calculating RG series for g 6 and applying Pade-Borel 
resummation technique to get proper numerical results. Actually, the first step in this 



direction has been already done by one of the authors [0 who has found two-loop RG 
expansion for g Q which yielded g$ = 1.50. This estimate, however, is undoubtedly very 
crude. Indeed, it has been obtained on the base of the lower-order perturbative expansion 
within the theory having no small parameter while accurate enough numerical estimates, as 
is well known, may be extracted only from sufficiently long RG series. Below we find the 
RG expansion for g§ in the four-loop approximation which will be shown to provide fair 
numerical estimate for the quantity of interest. 

The method of calculating g e we use here is straightforward. Since in three dimensions 
higher-order bare couplings are irrelevant in RG sense (see, e.g., Ref. 0]), renormalized 
perturbative series for g 6 can be obtained from conventional Feynman graph expansion of this 
quantity in terms of the only bare coupling constant - A. In its turn, A may be expressed 
perturbatively as a function of renormalized dimensionless quartic coupling constant g^. 
Substituting corresponding power series for A into original expansion we can obtain the RG 



series for g^. First two terms of this series have been presented earlier ||12|| . Thus, what we 
should find is the three-loop and four-loop contributions. As can be shown, they are formed 
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by 16 and 94 one-particle irreducible Feynman graphs, respectively. Their calculation gives: 
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The perturbative expansion for A emerges directly from the normalizing condition A = 
mZ 4 Z _2 g 4 and the well-known expansion for Z 4 : 
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Z 4 = 1 + — g 4 + —g\ + 1.778667825^ + O(gt). (4) 

Combining these expressions we obtain 

9e = -gl(l ~ -94 + 1.38996295104 - 2.50173240c/ 3 ). (5) 

Being a field-theoretical perturbative expansion, this series is divergent (asymptotic). 
Hence, direct substitution of the fixed point value g\ into fl5|) would not lead to satisfactory 
results. To get reasonable numerical estimate for g§ some procedure making this expansion 
convergent should be applied. Since the series ([5]) is alternating simple Pade-Borel technique 
may play a role of such a procedure. Pade approximants of [L/l] type, when used for 
analytical continuation of Borel transforms, are known to provide rather good results for 
various Landau- Wilson models (see, e.g., Refs. P,|5] flO|JTl|] ) . With the expansion (|5|) in hand, 



we can construct Pade approximant [2/1] for its Borel transform F(y) which is related to 
the function to be found ("sum of series") by the formula 
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f{x) = £ c k x k = [ e-*F(xt)dt, F(y) = £ % k . (6) 

k=0 J Q k=0 K - 

Numerical estimate for g§ is then obtained by evaluation of the Borel integral under g 4 = g\ 
where high-precision fixed point value of g 4 is known from six-loop RG calculations: g\ = 
0.988 The final result is as follows: 

9*6 = 1-596. (7) 
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How close to the exact value of g$ may this number be? To clear up this point let's 
compare the estimates given by four subsequent RG approximations. One-, two-, three-, 
and four-loop calculations of g§ give 2.763, 1.500, 1.622, and 1.596, respectively (the three- 
loop estimate has been presented earlier in Ref. |T3]]). Since this set of numbers demonstrates 
oscillatory (and rapid!) convergence one may expect that exact sextic effective coupling 
constant lies between the three-loop and four-loop RG estimates. So, our four-loop RG 
analysis leading to the number (0) underestimates gl by less than 2%. Moreover, it is 
possible to further improve the estimate for g§ by addressing the average 

1.622 + 1.596 

^ = 1.609 « 1.61, (8) 

which can be referred to as a most accurate approximate value differing from the exact one 
by no more than 1%. 

Another source of information about the accuracy of numerical results is their sensitivity 
to the type of resummation procedure. To get such an information we calculate g§ using the 
Borel-Leroy transformation 

oo 

f(x) = J t 6 e -'F(xt)rft, (9) 



which contains a free parameter b. This parameter is chosen in a way that provides fastest 
possible convergence of the iteration scheme. Being applied to the RG expansion (Q) the 
machinery described gives g$ = 1.61 under the optimal value of b. This number is seen to 
coincide with the estimate (§) what may be considered as a strong argument in favour of its 
high accuracy. 

The value of g\ we have just determined turns out to differ substantially from those 
obtained in Refs. |],@,|3|]. At the same time, it agrees fairly well with the estimate which 
follows from the three-loop e-expansion for the ratio g%jg\ ||15|| : 

9_1 = 2e - —e 2 + 1.2759e 3 + 0(e 4 ). (10) 
94 27 



Indeed, resumming this expansion by the Pade-Borel method and then putting e = 1 and 
#4 = 0.988 we find g£ = 1.653, i.e. the number which is sufficiently close to (H). Good 
agreement takes place also between our results and very fresh estimates g$ = 1.63 ±0.05 and 
gl = 1.57 ± 0.10 deduced from strong coupling series and high temperature expansions for 
the lattice model [16]. At last, numbers (0), (§) are in a perfect accord with another very 
recent estimate g$ = 1.604 |17| found by the analysis of the five-loop scaling equation of state 
|T8| which has used more sophisticated than ours resummation technique (the Borel-Leroy 
transformation combined with a conformal mapping). 

To summarize, we have calculated the four-loop RG expansion for dimensionless sextic 
effective coupling constant g$ of the three-dimensional Ising model. Resummation of this 
expansion by the Pade-Borel method has lead at criticality to the result g$ = 1.596. Having 
analyzed the set of four subsequent RG estimates for g$ the number 1.61 has been argued 
to be the best approximate value of universal sextic effective coupling accurate within 1%. 
Precisely the same estimate has been found by applying another, Pade-Borel-Leroy proce- 
dure to resum the four-loop RG series. Such a stability of the obtained value under the 
change of resummation technique may be considered as a serious argument in favour of its 
high accuracy. Although g$ = 1.61 is substantially smaller than the early estimates |],[2||| 
it turns out to be in a good agreement with that following from the Pade-Borel resummed 



e-expansion and with the results of recent advanced calculations |T6|,|T7 
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